The Anjials of Probability 
2011, Vol. 39, No. 5, 1844-1863 
DOI: 10.1214/11-AOP652 

@ Institute of Mathematical Statistics, 2011 

SCHRAMM'S PROOF OF WATTS' FORMULA 

By Scott Sheffield and David B. Wilson 
Massachussetts Institute of Technology and Microsoft Research 

Gerard Watts predicted a formula for the probability in percola- 
tion that there is both a left -right and an up-down crossing, which 
was later proved by Julien Dubedat. Here we present a simpler proof 
due to Oded Schramm, which builds on Cardy's formula in a concep- 
tually appealing way: the triple derivative of Cardy's formula is the 
sum of two multi-arm densities. The relative sizes of the two terms 
are computed with Girsanov conditioning. The triple integral of one 
of the terms is equivalent to Watts' formula. For the relevant calcula- 
tions, we present and annotate Schramm's original (and remarkably 
elegant) Mathematica code. 

1. Watts' formula. When Langlands, Pichet, Pouliot and Saint-Aubin 
(1992) were doing computer simulations to test the conformal invariance 
of percolation, there were several different events whose probability they 
measured. The first event that they studied was the probability that there 
is a percolation crossing connecting two disjoint boundary segments. Us- 
ing conformal field theory, Cardy (1992) derived his now-famous formula 
for this crossing probability, and the formula was later proved rigorously 
by Smirnov (2001) for site percolation on the hexagonal lattice. The next 
event that Langlands et al. tested was the probability that there is both 
a percolation crossing connecting the two boundary segments and a per- 
colation crossing connecting the complementary boundary segments (see 
Figure 1). This probability also appeared to be conformally invariant, but 
finding a formula for it was harder, and it was not until several years after 
Cardy's formula that Watts (1996) proposed his formula for the probabil- 
ity of this double crossing. Watts considered the derivation of the formula 
unsatisfactory, even by the standards of physics, but it matched the data 
of Langlands et al. very well, which lent credibility to the formula. Watts' 
formula was proved rigorously by Dubedat (2006a). 
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Fig. 1. In the left panel, there is no left-right crossing in blue hexagons. In the second 
panel there is a blue left-right crossing, but no blue up-down crossing. In the third panel, 
there are both blue left-nght and blue up-down crossings. Cardy's formula gives the prob- 
ability of a left-right crossing in a domain, while Watts ' formula gives the probability that 
there is both a left-right crossing and an up-down crossing. 



To express Cardy's formula and Watts' formula for the two types of cross- 
ing events, since the scaling limit of percolation is conformally invariant, it 
is enough to give these probabilities for one canonical domain, and this is 
usually taken to be the upper half-plane. There are four points on the bound- 
ary of the domain (the real line). Label them in increasing order xi, X2, X3 
and 2:4. Cardy's formula is then the probability that there is a percolation 
crossing from the interval [xi,X2] to the interval [x3,X4]. Again by conformal 
invariance, we may map the upper half-plane to itself so that xi — )■ 0, X3 — )■ 1 
and X4 — >• 00. The remaining point X2 gets mapped to 

n IN f \ {X2-Xi){x4-X3) 

(1.1) s = cr{xi,X2,xs,X4) := -, 

(Xs -Xi)[X4 -X2) 

which is a point in (0, 1) known as the cross-ratio. Both Cardy's formula and 
Watts' formula are expressed in terms of the cross-ratio. Cardy's formula for 
the probability of a percolation crossing is 

r(2/3) „i/3,_,.Yi 2.4, 



(1.2) cardy(s) := , , \ / , , s^^'^ 2F1 -,-;-; s 

^ ' ' r(4/3)r(l/3) V3 3 3' 

where T is the gamma function, and 2F1 is the hypergeometric function 
defined by 

77/1, \ {o)n{V)n n 

n.=0 

where a, 6, c € C are parameters, c ^ — N (where N = {0, 1,2,.. .}), and (i)n 
denotes i{i + 1) • • • (^ + n — 1). This series converges for z G C when \z\ < 1, 
and the hypergeometric function is defined by analytic continuation else- 
where (though it is then not always single- valued). 

By comparison, Watts' formula for the probability of the two crossings is 
the same as Cardy's formula minus another term 

r(2/3) ^1/3 ^fl 2_4_ 



watts(s) := — — , — —s ' o-^i I — , — ; —',s 
^ ' r(4/3)r(l/3) V3 3'3' 
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FlG. 2. Cardy's formula (upper curve), Watts' formula (lower curve), and a tripod prob- 
ability (defined in Section 3) as a function of the cross-ratio s. 



where 3F2 is the generahzed hypergeometric function. The functions cardy(s) 
and watts(s) are shown in Figure 2. [The reader should not be intimidated 
by these formulae; the parts of the proof involving hypergeometric func- 
tions can be handled mechanically with the aid of Mathematica. See also 
Watts (1996) and Maier (2003) for equivalent double-integral formulations 
of Watts' formula.] 

Schramm thought that Dubedat's paper on Watts' formula was an excit- 
ing development and started reading it as soon as it appeared in the arXiv. 
Schramm sometimes presented papers to interested people at Microsoft Re- 
search: for example, he presented Smirnov's proof of Cardy's formula when 
it came out [Smirnov (2001)], as well as Dubedat's paper on Watts' for- 
mula [Dubedat (2006a)], and later Zhan's paper on the reversibility of SLE^ 
for K < 4 [Zhan (2008)]. In the course of reaching his own understanding 
of Watts' formula, Schramm simplified Dubedat's proof, with the help of 
a Mathematica notebook, and it was this version that he presented at Mi- 
crosoft on May 17, 2004. This proof did not come up again until an Au- 
gust 2008 Centre de Recherches Mathematiques (CRM) meeting on SLE 
in Montreal, after a talk by Jacob Simmons on his work with Kleban and 
Ziff on "Watts' formula and logarithmic conformal field theory" [Simmons, 
Kleban and Ziff (2007)]. Schramm mentioned that he had an easier proof 
of Watts' formula, which he recalled after just a few minutes. The people 
who saw his version of the proof thought it was very elegant and strongly 



(1.3) 
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encouraged him to write it up. The next day Oded wrote down an outhne 
of the proof, but he tragically died a few weeks later. There is interest in 
seeing a written version of Schramm's version of the proof, so we present it 
here. 

2. Outline of proof. This is a slightly edited version of the proof outline 
that Oded wrote down at the CRM. Steps 1 and 2 are the same as in 
Dubedat's proof, but with step 3 the proofs diverge. We will expand on 
these steps of the outline (with slightly modified notation) in subsequent 
sections. 

• Reduce to the problem of calculating the probability that there is a cross- 
ing up-down which also connects to the right. 

• Further reduce to the following problem. In the upper half-plane, say, mark 
points — oo < yi < xo < 2/2 < 2/3 = oo. Let 7 be the SLEg interface started 
from xq. Let r := mi{t > : 74 € M \ [yi, 1/2] and a := sup{i < r : 7^ G M}. 
Calculate P[7r € [2/2, 2/3], Ta G [2:0, 2/2]]- 

• Let fJi := sup{t <T:^t ^ [yi^^o]} and (72 := sup{t KT-.jt^ [xo,y2]}- Now 
calculate the probability density of the event 70-1 = zi , 70-3 = ^2 , 7t = -23 as 
h{zi,xo,Z2,Z3) := dz^dz^dzs Cardy(zi, xq, 22, ^3)- 

• Now, h [times certain derivatives] is a martingale for the corresponding 
diffusion. Consider the Doob-transform (/i-transform) of the diffusion with 
this h. This corresponds to conditioning on this probability zero event. For 
the Doob-transform, calculate the probability that (T2 > ai. This comes 
out to be a hypergeometric function g. Finally, 



(or more precisely, the three-arm probability), and use integration by 



3. Reduction to tripod probabilities. The initial reduction, which is step 1 
of the proof, has been derived by multiple people independently. The first 
place that it appeared in print appears to be in Dubedat's (2004) paper, 
where it is credited to Werner, who, in turn, is sure that it must have been 
known earlier. In the interest of keeping the exposition self-contained, we 
explain this reduction. 

It is an elementary fact that exactly one of the following two events oc- 
curs: 

(1) there is a horizontal blue crossing in the rectangle (i.e., a path of blue 
hexagons connecting the left and right edges of the rectangle), which we 
denote by Hi,; 

(2) there is a vertical yellow crossing, which we denote by Vy. 




parts. 
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If there is a horizontal crossing, then by considering the region beneath 
it, using the above fact, either it connects to the bottom edge (forming a T 
shape) or else there is another crossing beneath it of the opposite color. Since 
there are finitely many hexagons, there must be a bottom-most crossing, 
which then necessarily forms a T shape. Thus exactly one of the following 
three events occurs: 



(1) there is no horizontal crossing of either color (denoted by A^); 

(2) there is a blue T (denoted Th); 

(3) there is a yellow T (denoted Ty). 




Of course the latter two events have equal probability, so we have 

Pr[iV] +2Pr[rfe] = 1. 

Recall again that there is either a blue horizontal crossing or a yellow ver- 
tical crossing but not both. We can decompose the yellow vertical crossing 
event into two subevents according to whether or not there is also a yel- 
low horizontal crossing. The first subevent is, of course, the event we are 
interested in (with blue and yellow reversed), and the second subevent is 
identical to the event A'^. 




Thus we have 



Pv[Hb] + PT[Hy A Vy] + Fi[N] = 1. 

Combining these equations, we see that 

PT[HbAVb]=2Fr[n]-Fr[Hb]. 

In the limit of large grids with cross ratio s, the third term is given by 
Cardy's formula, cardy(s), and we seek to show that the left-hand side is 
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given by Watts' formula, watts(s). Let us give another name for what we 
expect to be the limit of the second term. Define tripod(s) to satisfy 

watts(s) = 2tripod(s) — cardy(s), 

that is [substituting (1.2) and (1.3)], 

watts(s) + cardy(s) 



tripod(s) 



2 

r(2/3) 



r(4/3)r(i/3) 



1 / 4 5 

-S3F2 l,l,-;2,-;s 



2r(l/3)r(2/3) V ' '3' '3 

Then in light of Cardy's formula, proving Watts' formula is equivalent to 
showing that Pr[T(,] is given by tripod(s) in the fine mesh limit. 



4. Comparison with SLEe. 



4.1. Discrete derivatives of the tripod probability. Consider percolation 
on the upper half-plane triangular lattice, and let 2:2, X3, 2:4] be the 

probability of a blue tripod connecting the intervals (xi, ^2) and (x2, 2:3) and 
(2^3, 2:4) when the four (here discrete) locations are xi < X2 < X3 < X4 (each 
of which is a point between two boundary hexagons; see the upper image in 
Figure 3). 

Then A^;^ Pt[xi,X2,X3,X4\ := PT[xi,X2,X3,Xi\- Pt[xi,X2,X3,X4-1] gives 
the probability that there is a crossing tripod for (xi, 0:2, 2:3, X4) but not for 
{xi,X2,X3,X4 — 1). (Here we assume that the lattice spacing is 1.) Since 
the crossing tripod for (xi, 3:2, X3, 2:4) does not extend to a crossing tripod 
for {xi,X2,X3,Xi — 1), there must be a path of the opposite color from the 
hexagon just to the left of X4 — 1 to the interval between a;2 + 1 and X3 + 1; 
this event is represented by the second image in Figure 3. Similarly, 

-A^-^A^^A^^PT[xi,X2,X3,Xi\ 

gives the probability of a multi-arm event such as the one in the bottom 
image in Figure 3. 

By summing these discrete differences, it is straightforward to write 

Pt[xi,X2,X3,X4\= X] X] -A^^Ax.^A^^PT[a,X2,b,c\. 

c£{x-i,X4] b£{x2,x-i] a£{xi,X2] 

If there is a blue tripod connecting the intervals {xi,X2), ix2,X3) and (x3,X4), 
then there is only one cluster containing such a tripod. This formula can be 
interpreted as partitioning the tripod event into multi-arm events of the 



SCHRAMM'S PROOF OF WATTS' FORMULA 



7 




Fig. 3. The discrete triple partial derivative of the tnpod probability is the probability of 
a multi-arm event. The top panel illustrates the event whose probability is Pt[xi,X2,X3,X4], 
the next panel illustrates Ax^Pt[xi,X2,X3,X4\, the third panel illustrates Ax^Ax4Pt[xi,X2, 
X3,X4,] and the bottom panel illustrates —AxiAx3Ax^Pt[xi,X2,X3,X4]. 

type shown in the bottom panel of Figure 3. The triple {a,b,c) G {xi,X2] x 
(x2,X3] X {x3,X4] is uniquely determined by the tripod: a is (half a lattice 
spacing to the right of) the rightmost boundary point of the tripod cluster 
in the interval {xi,X2), b is (just right of) the rightmost point of the tripod 
cluster in (2:2, X3) and c is (just right of) the leftmost point of the tripod 
cluster in (X3, X4). 

4.2. Discrete derivatives of the crossing probability. Consider percola- 
tion on a half-plane triangular lattice, as in the previous subsection, and 
let Pc[xi,X2,X3,X4] be the probability of at least one blue cluster span- 
ning the intervals {xi,X2) and (2^3, X4); see the upper image in Figure 4. 
Then A^^Pc[xi,X2,X3,X4] = Pc[xi,X2,X3,X4] - Pc[xi,X2,X3,Xi - 1] gives 
the probability that there is a crossing for {xi,X2,X3,X4) but not {xi,X2,X3, 
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Fig. 4. The discrete triple partial derivative of the crossing probability is sum of the 
probabilities of two multi-arm events. The panels illustrate the events whose probability 

is Pc[xi,X2,X3,Xi] (top), Aa:^Pc[Xl,X2,X3,X4] (sCCOud), Aj:^Aj:_^Pc[xi,X2,X3,X4\ (third 

row), and —A.j:-^^Ax3Ax^Pc[xi,X2,X3,Xii\ (bottom row). 

X4 — 1). This event is represented by the second image in Figure 4. Similarly, 

- Ax-i Aa;3 A Pc , 2:2 , , 3:4] 

gives the probability that one of the two multi-arm events in the bottom 
image in Figure 3 occurs. The event of a crossing cluster is equivalent to the 
event that one of these multi-arm events occurs for some (necessarily unique) 
set of three points (a, 6, c) € (xi, 3:2] x (3:2, X3] x (xs, X4]: a is (just right of) the 
rightmost boundary point of the crossing cluster (s) in the interval (xi,X2); 
b is (just right of) the rightmost point of the crossing cluster in (x2,X3) [if it 
exists; otherwise h is the rightmost boundary point in (x2,X3) of a crossing 
yellow cluster, as shown]; and c is (just right of) the leftmost point of the 
cluster (s) in (x3,X4). 

Thus — Aa;^ Aa;3 Aa;^Pc'[xi, X2, X3, X4] decomposcs into the probabilities of 
two multiarm events, the first of which is — A^;^ Aa;3Aa;4Pr[xi, X2, X3, X4]. 
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Fig. 5. The interface interpretation of the multi-arm events. 



4.3. Multi-arm events and the interface. Consider the setting of Fig- 
ures 3 and 4, and suppose we add an additional boundary layer of blue 
hexagons to the left of X2 and yellow hexagons to the right of X2- Then 
let 

Tdiscrete be the discrete interface starting at X2- (See Figure 5.) 
Then the union of the two multi-arm events at the bottom of Figure 4 de- 
scribes the event that that c is the first boundary point that 7discrcte hits 
outside the interval (xijX^), and that a and b are the leftmost and rightmost 
boundary points hit by 7discrctc before c. The left figure corresponds to the ca- 
se that a is hit before 6, and the right figure to the case that b is hit before a. 

4.4. Continuum Watts ' formula: A statement about SLE. Like Cardy's 
formula, Watts' formula has a continuum analog, which is a statement 
strictly about SLEg. Fix real numbers xi < X2 < X3 < X4, and let s be their 
cross ratio. Consider the usual SLEg in the upper half-plane, where the 
starting point of the path is X2- Before Smirnov proved Cardy's formula for 
the scaling limit of triangular lattice percolation, it was already known by 
Schramm that cardy(s) represents the probability that 7 hits [x^^X/^) before 
hitting R \ [ xi,X4]. [In the discrete setting of Section 4.3, having 7discrete 

hit 

(x3,X4) before the complement of (xi,X4) is equivalent to the existence of 
a crossing.] In light of Section 4.3, the following is the natural continuum 
analog of the tripod formula. 

Theorem 4.1. Let SLEtripod(s) be the probability that both: 

(1) 7 first hits {x^^Xij (at some time t) before it first hits ]R\ {xi,X4^), 
and 

(2) 7 hits the leftmost point o/Mn7[0,t) before it hits the rightmost point. 
Then SLEtripod(s) = tripod(s). 
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Theorem 4.1 is the actual statement that was proved by Dubedat, and the 
statement whose proof was sketched by Oded. Dubedat claimed further that 
Theorem 4.1 would imply the tripod formula (and hence Watts' formula) for 
the scaling limit of critical triangular lattice percolation if one used the (at 
the time unpublished) proof that SLEg is the scaling limit of the interface 
[Dubedat (2006a)]. To be fully precise, one needs slightly more than the fact 
that the interface scaling limit is SLEg: it is important to know that the 
discrete interface is unlikely to get close to the boundary without hitting it. 
[Similar issues arise when using Cardy's formula to prove SLEg convergence; 
see, e.g., Camia and Newman (2007).] Rather than address this (relatively 
minor technical) point here, we will proceed to prove Theorem 4.1 in the 
manner outlined by Oded and defer this issue until Section 7. 

It is convenient to have a name for the SLE versions of the multi-arm 
events in Figure 4. Say that a triple of distinct real numbers (a, 6, c) with 
a < < 6 constitutes a tripod set for 7 if for some t> we have: 

(1) 7(i) = c; 

(2) inf(7[0,t) nM) =a; 

(3) sup(7[0,t)nM) = 6. 

There are a.s. a countably infinite number of tripod sets, but if xi < and 
X3 > is fixed, there is a.s. exactly one for which xi<a<0<6<X3 and 
c ^ (xi,X3). Let Vc'^^ — )■ M be the probability density function for this 
{a,b,c). (We see in Lemma 4.2 that this density function exists.) There are 
also two types of tripod sets (a, 6, c): those for which 7 hits a first and those 
for which 7 hits b first. Write Vc = T^A + "Pb, where Va and Vb are the 
corresponding probability densities for a- first and 6- first tripod sets. 
Now, we claim the following: 

Lemma 4.2. Using the notation above, the density functions Vc and Va 
exist, and 

cardy(s) = / / / Vc{a,b,c) dcdbda 

J xi Jo J X3 

and 

rO rxs f-Xi 

SLEtripod(s) = III VA{o,,b,c) dcdbda. 

J xi Jo J X3 

Proof. In the event that the density functions do not exist, we abuse 
notation and let ^^Vc{a,b,c) dcdbda^^ and ^^VAia,b,c) dcdbda^^ denote the 
relevant measures, which must exist. It is easy to see that the event that 
{x3,X4) is hit before M\ [2;i,X4] is equivalent to the event that (a,6, c) G 
(xi,0) X (0,^3) X (2:3, X4). By definition, SLEtripod(s) is the probability 
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of the same event intersected with the event that 7 hits a first. Finally, 
observe that Cardy's formula of the cross-ratio of 4-points is three-times 
differentiable, so the density function Vc{a-,h,c) exists and, consequently, 
the density function VA{0'^b^c) also exists. □ 

Of course, from this, one has the immediate corollary: 

Corollary 4.3. Using the above notation, 

dx^dx^dx^ cardy(cr(xi, 0, X3, X4)) = Vcixi^x^, X4) 

and 

dxi dxa 0x4, SLEtripod(cr (xi , 0, X3 , X4) ) = Pa (xi , X3 , X4 ) . 

If we could show further that 

(4.1) 'Pa{xi,X3,X4,) = 8x^8x38x4 tripod(cr (xi , 0, X3 , X4) ) , 

then this corollary and standard integration would imply Theorem 4.1, 
since we know that tripod(cr(-)) = SLEtripod(cr(-)) on the bounding planes 
xi = 0, X2 = and X3 = X4. Since we already have an explicit formula for 
tripod, the only remaining step is to explicitly compute Va- Oded's ap- 
proach is to compute the ratio Va/'Pc as the conditional probability [given 
that {a,b,c) form a tripod set] that 7 hits a before b. Since Vc is known, 
this determines Va- 

5. Conditional probability that a is hit first. Schramm was very adept 
with using Mathematica to calculate all manner of things. He probably 
would have considered this last step to be routine, since it was for him 
straightforward to set up the right equations and then let Mathematica 
solve them. At this point we refer to his original Mathematica notebook 
from 2004, and explain the various steps in the calculation. To be consistent 
with Oded's notation, we now make the following substitutions: 

V3 = a, W = X2, vi = b, V2 = c. 

(We assume V3 <W < vi < V2. Oded apparently chose this notation because 
under cyclic reordering it was the same as W = vo,vi,V2,V3.) 

First we formally define the function cardy(s) as in (1.2). The Mathe- 
matica function cardy defined here involves an additional parameter k, but 
it specializes to cardy(s) when k = 6. This more general formula is analo- 
gous to Cardy's formula but gives the (conjectural) crossing probability for 
the critical Fortuin-Kasteleyn random cluster model [with q = 4cos^(47r/K)] 
with alternating wired-free-wired-free boundary conditions. (See Rohde and 
Schramm [(2005), conjecture 9.7], for some background.) This formula was 
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known at Microsoft in 2003, and most likely Oded copied it from another 
Mathematica notebook. This formula was later independently discovered 
[Bauer, Bernard and Kytola (2005)] (nonrigorously) and [Dubedat (2006b)] 
(rigorously) . 



l_i r-4+K 4 -4+K 1 

oardy[s_] = c [x] s x Hypergeoinetric2Fl , — , 1 + , s 



l_i r-4+K 4 -4+K 1 

s K c [k] HYpergeometric2Fl , — , 1 + , s 

K K K ^ 



Consider the evolution of chordal SLEg started from W and run to oo, 
when at time zero there are 3 marked points at positions vi, V2 and ^3. We 
then let W{t) represent the SLEg driving function [i.e., W{t) = W{Qi) + \/^Bt 
where Bt is a standard Brownian motion] of the Loewner evolution 

and interpret the Vi as functions of t, evolving under the Loewner flow, that 
is, Vi{t) := gt{vi{0)). 

If / is any function vi,V2,V3,W , we define 

(5.1) L{f) := ^E[f{W{t)Mt),V2{t)Mm 

This is a new function of the same four variables which can be calculated 
explicitly using Ito's formula as 

^2 3 



1=1 

This operator is defined as L in the Mathematica code below. 



L[f_] := D[f, W, W] K/2 + Suin[D[f. v[i] ] 2 / (v[i] - W) . {i, 1, 3}] 



^ (W-v[3]) (v[l] -v[2]) 
(W-v[2]) (v[l] -v[3]) 

(v[l] -v[2] ) (W- v[3] ) 
(W-v[2]) (v[l] -v[3]) 



Similarly in the Mathematica code, cr is the cross-ratio [defined in (5.2)], 
that is, 

(5.2) cT:=cr{v3,W,vi,V2) = --r-. r. 

[W - V2){Vi - V3) 
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In the next line, Oded performed a consistency check. Cardy's formula 
should be a martingale for the SLE^ diffusion, hence L(cardy(cr(-))) = 0. 



L[cardy[cr]] /. k -» 6 // Simplify 



Next, Oded computes the triple derivative h of Cardy's formula (which is 
the same as the Vc defined in Section 4.4). 

h = D[cardy[cr] , v[l] , v[2] , v[3] ] /. k ^ 6 // Simplify 

(2 c [6] (W - v[l] ) (W - v[3] ) 

(v[2]2v[3]2 -v[l] v[2] v[3] (v[2] +v[3]) + 
v[l]2 (v[2]2 - v[2] v[3] + v[3]2) + 

(v[l]2+v[2]2-v[2]v[3] + v[3]2-v[l] (v[2] + v[3])) - 
w(v[l]2 (v[2] + v[3]) + v[2]v[3] (v[2] + v[3]) + 
v[l] (v[2]2-6v[2] v[3] +v[3]2))))/ 
(27 (W-v[2])'' (v[l] -v[3])5 (((W-v[l]) (v[l] -v[2]) 

(W-v[3]) (v[2] -v[3] ) ) / ( (W- v[2] )2 (v[l] - v [ 3 ] ) ^ ) ) ) 



That is, he computes 



The result is somewhat complicated, but we may ignore it, since it is just 
an intermediate result. 

The next step involves conditioning on an event of zero probability, the 
event that {v^^vi,V2) is a tripod set. We can make sense of this by intro- 
ducing a triple difference of Cardy's formula and recalling the results of 
Section 4.4. First we introduce notation to describe some small evolving 
intervals. For given values ui(0), 1)2(0), f 3(0), VF(0), pick e small enough so 
that the intervals (iJj (0) , (0) + e) are disjoint and do not contain 1/F(0). 
Write Vi{Q) = Vi(0) + e. Define Vi{t) using the Loewner evolution, and write 
ei{t) ■.= Vi{t) —Vi{t). Let us write 









where A^, is the difference operator defined by 



/:^Vf{v) = f{v + e)-f{v). 
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Note that the A^,'^'^ depend on t. By Corohary 4.3, equation (5.3) at thne t 
represents the conditional probabihty (given the Loewner evolution up to 
time t) that there is a tripod set in [v^ (0) , -ua (0)] x [vi (0) , ui (0)] x [^;2 (0) , 'i)2 (0)] . 
By Girsanov's theorem, conditioning on this event induces a drift on the 
Brownian motion Wt driving the SLE, where the drift is 

K dw log ,£2,£3 (^^3 ,W,Vi,V2). 

Observe that 

_d_ _ {d/dw)h,,,,,,,, _ !jj{d/dw)h 

5^iOg/.,,,,„,3- ^^^^^^^^^ - ^^^^ 

~ r = dw log h, 

£ie2e^h 

where there triple integral is over j, -Sj]. Thus upon taking the limit e — )• 0, 
the drift becomes 

(5.4) drift(t) := k5vi/ log /i(u3, VF, fi, f2)- 

The next Mathematica code explicitly computes (5.4). 

drift = xD[Log[h] , W] /. 6 // Simplify 

(2 (-2v[2]^v[3]^ + 3v[l]v[2]2v[3]2 (v[2] + v[3]) + 
3v[l]2v[2]v[3] (v[2]2-4v[2]v[3] + v[3]2) + 
v[l]^ (-2v[2]^ + 3v[2]2v[3] + 3v[2]v[3]2-2v[3]^) + 

(-2v[l]^-2v[2]3 + 3v[2]2v[3] + 3v[2]v[3]2-2v[3]^ + 
3v[l]^ (v[2] + v[3]) + 3v[l] (v[2]2-4v[2]v[3] + v[3]2)) + 
3W2 (v[l]^ (v[2] + v[3]) + v[l]2 (-4v[2]2 + 2v[2] v[3] - 
4v[3]2) + v[2] v[3] (v[2]^-4v[2] v[3] + v[3]2) + 
v[l] (v[2]3 + 2v[2]2v[3] + 2v[2]v[3]2+v[3]^)) + 
3w(v[2]2v[3]2 (v[2] + v[3]) + 2v[l]v[2]v[3] 
(-2v[2]2+v[2]v[3] - 2v[3]2) + 
v[l]^ (v[2]2-4v[2]v[3] + v[3]2) + 

v[l]2 (v[2]3 + 2v[2]2v[3] +2v[2] v[3]2+v[3]3))))/ 
((W-v[l]) (W-v[2]) (W-v[3]) (v[2]2v[3]2- 
v[l] v[2] v[3] (v[2] + v[3] ) + 
v[l]2 (v[2]2 -v[2] v[3] +v[3]2) + 

(v[l]2+v[2]2-v[2]v[3] + v[3]2-v[l] (v[2] + v[3])) - 
W (v[l]2 (v[2] + v[3]) + v[2] v[3] (v[2] + v[3]) + 
v[l] (v[2]2-6v[2] v[3] +v[3]2)))) 
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The expression above is complicated, but again it is an intermediate result 
that we do not need to calculate or read ourselves. The first line of the 
Mathematica code below defines the generator LI (which we will write as Li) 
for the conditioned SLEg, where the driving function Wt has the drift given 
above. Here Li is defined as in (5.1) except that the expectation is with 
respect to the law of Wt with the drift term (5.4). Thus 

Li{f):=L{f)+dTimdwf. 

As before, if / is a real function of W, vi,V2,V3, then Li{f) will be a function 
of the same four variables. 

We now compute the probability in the modified diffusion that is 
absorbed before vi, that is, that W{t) collides with vs{t) before colliding 
with vi (t) . This probability will be a martingale that only depends upon the 
cross-ratio s. Thus, in the next paragraph, we specialize and consider func- 
tions of W,vi,V2,V3 that have the form f{cr{v3,W,vi,V2)) where /:R— )-R 
is a function of one variable. We would like to find a one-parameter func- 
tion / for which f{cr{v3, W, 1^1,^2)) is a martingale with respect to this mod- 
ified diffusion, so we will require that Li{f{cT{v3,W, vi,V2))) = 0. What one- 
parameter functions / have this property? 

Oded answers this question with some clever Mathematica work. First, he 
re-expresses the differential equation Li{f{cr{v3,W,vi,V2))) = — which in- 
volves the four parameters W,vi,V2,vs — in terms of the parameters s,vi,V2, 
V3. He does this by setting s equal to the expression for cr given in (5.2), 
solving to get W in terms of the other variables and plugging this new 
expression for W into the expression Li(/(cr(r;3, VK, ui, U2))). 



Ll[f_] := L[f] +driftD[f, W] 



Ll[f [cr] ] /. Solve [cr == s, W] [ [1] ] /. k 6 // Simplify 

(((-1 + s) v[l] +v[2] -sv[3])'' 

(2 (1-6 + 4 s^) f'[s] + 3s (-l + 2s-2s^ + s^) f"[s]))/ 
(s (-l + 2s-2s2+s3) (v[l] -v[2])2 (v[l] -v[3])2 (v[2] -v[3])2) 



This expression for Li{f{cic{v3,W,vi,V2))) depends on /' and /", and 
equating it to zero yields a differential equation for /, the unknown one- 
parameter function of the cross-ratio that we seek, 

2(1 - 6s2 + 4s3)/'(s) + 3s(-l + 2s- 2s^ + s^)f"{s) = 0. 

Oded solves this differential equation, which yields the function / up to two 
free parameters Ci and C2. 
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DSolve [% == 0, f [s] , s] 


{{f [s] ^ C[2] + 

(1 - s 

(3 (-1 + s) 


's^--^ c[l] 
(1-s 

(l_s + 


l-3s + 2s2- 

f ) HYpergeometric2Fl ^ — , — , — , sjjjy^ 

=•))}} 



Here Mathematica gives 

,^3 i-3s + 2s^-{l- -s + s^) 2Fi(2/3, 1/3; 5/3; s) 

3(— 1 + s)^'-^[l — s + s'') 

The conditional probability that we seek tends to 1 when s — )■ and tends 
to when s — t- 1, and this determines Ci and C2: C2 must be 1, and Ci 
follows from Gauss's hypergeometric formula, 

, , r(c)r(c-a-b) 

2F1 a,b;c;l = ' \ 

r(c — a)T[c — b) 

Solving for Ci and substituting, we find that the conditional probabihty 
that 7 hits V3 before vi, given that {v3,vi,V2) is a tripod set, is given by 



r(2/3)r(5/3) 



1 + 3s -2^2 „ /2 1 5 



(l-s)i/3(i_s + s2) - \3'3'3 
where s is the cross-ratio of vs,0,vi,V2- 

6. Comparison of triple derivatives. Taking Va and Vc as defined in 
Section 4.4, and / and h as defined in the previous section, we now have 

h{v3,0,Vi,V2) =Vc{v3,Vl,V2) 

and 

f{cT:{v3,0,Vi,V2))=VA{v3,Vi,V2)/Vc{v3,Vl,V2). 

In principle the next step toward proving (4.1) (and hence Theorem 4.1) 
would be to integrate Va = /(cr(-))/i over the three variables 1^1,^2, "^3 and 
show that one obtains tripod(cr(-)). In Oded's original notes, he stated that 
this could be done using integration by parts. Fortunately (for those who 
lack Oded's skill at integrating) we already know (thanks to Watts) what 
we expect tripod to be, so we can instead differentiate tripod(cr(-)) three 
times (w.r.t. vi,V2,V3), and check that it equals /(cr(-))/i. The Mathematica 
code in this final section was generated by the authors of this paper, not by 
Schramm. 

First we redefine cardy to have an explicit constant and define the pur- 
ported tripod probability. 
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Gamma - 

caray[s_] = s^'^ 

Gamnia - J Gamma - J 


rl 2 4 1 

Hypergeometric2Fl — , — , — , s ; 

La 3 3 J 


tripod [s_] = 

cardy[s] - HypergeometricPFQ 1, — |, |2, — j, sj ; 

2 Gamma - j Gamma - j ^ ^ 


Next we differentiate Cardy's 


formula three times. 


dddc = D[cardy [cr] , v[l] , v[2] , v[3] 


] /. Solve [or s, W] [ [1] ] // FullSimplify 


2 22'3 


(1 + (-1 + s) s) 


27 (- (-1 + s) s)^'^ Gaimna[ij Gamma ^ ^■ 


(v[l] - v[2] ) (v[l] -v[3] ) (v[2] - v[3] ) 



This is the same as h defined earlier, but with the trick of eliminating 
the variable W and expressing the formula in terms of s. Next we triply 
differentiate the purported tripod probability. 



dddt = D[tripod[cr] , v[l] , v[2] , v[3] ] /. Solve[or - s. 


W] [[1] 


] // FullSimplify 


^4 (1+ (-1 + S) s) Gamma [ - ] + 






?n r 1 W r 1 2 5 
s Gamma — - ( 1 + s) HYpergeometrxc2Fl — , — , — , 
L 3 J I L 3 3 3 






r 2 4 5 

(-1 + s) (-l + 2s} Hypergeometric2Fl — , — , — , s 

L 3 3 3 






|3 6 sfY n (-(-1 + s) s)^'^ Gamma[^j (v[l] - v[2]) (v[l] - 


v[3] ) 


(v[2] -v[3])| 



Notice that the triple derivative of the tripod probability is expressed in 
terms of two different hyper geometric functions. In order to compare this ex- 
pression with the conditional probability computed in Section 5, we need to 
use some hypergeometric identities. We use one of Gauss's relations between 
"contiguous" hypergeometric functions [Erdelyi et al. (1953), Section 2.8, 
equation 33], to write 

-3 2-^i(i' h h ■5) + 1(1 - ■5) 2^i(f ' |; |; ^) - 5 2^i(i' f; ^) = 0- 

But 2Fi{c,b;c;s) = (1 — s)~^ Erdelyi et al. (1953), Section 2.8, equation 4, 
so 



(6.1) 



2Fi(|,|;|;.)=2(l-.)-V3_,F^(2^i;|;,). 
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dddt = 

dddt /. Hypergeometrio2Fl[2 / 3, 4/3, 5/3, s] -» 

2 (1 - s) * (-1 / 3) - Hypergeoinetrio2Fl [1 / 3, 2/3, 5/3, s] // Simplify 

ii (1+ (-1 + S) s) Gamma [ - 1 2 s^'^ Gamma [ ^1 

f-l + 3s-2s^+ (l-s + s^) HYpergeometric2Fl \-, - , - , si 1 

\ ^ ' L 3 3 3 J J _ 

iie -■JY 7T (-(-l + s) s)2" Gamma[-j (v[l] - v[2]) (v[l] - v[3]) (v[2] - v[3]) 



Next we compare the two expressions for the conditional probabihty and 
verify that they are the same. 



CP = 


dddt/dddc; 




cp2 








1 




Gamma - 

Gamma ^ - j Gamma ^ - j 

-l + 3s-2s2 r215i 
+ Hypergeometric2Fl — , — , — , s 

{l-s)i/3 (i_s + s2) L3 3 3 J 




cp - 


cp2 // FullSimplify 








Therefore the triple derivatives agree, and we have established (4.1). 



7. Percolation statement. We have the established equivalence of 
SLEtripod(s) and tripod(s), but we still need to make the connection to 
percolation. 

Theorem 7.1. Let D cC be a fixed bounded Jordan domain with marked 
points xi,X2,xs,X4 on its boundary. For any e, we may consider the hexag- 
onal lattice rescaled to have side length e and color the faces blue and yellow 
according to site percolation. Let B he the closure of the set of blue faces, 
and let be the probability that BOD contains a connected component that 
intersects all four boundary segments (xi,X2), (x2,X3), (x3,X4) and (x4,xi). 
Then 

lim P*^ = watts (s). 

Proving Theorem 7.1 solves the problem addressed by Watts. However, we 
remark that more general statements are probably possible. Any domain D 
with four marked boundary points has a "center" c{D) with the property 
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that a conformal map taking the domain to a rectangle (and the points 
to the corners) sends c to the center of the rectangle. Oded would proba- 
bly have preferred to show that for any sequence D„ of simply connected 
marked hexagonal domains (domains comprised of unions of hexagons within 
a fixed hexagonal lattice H with four marked boundary points of cross ra- 
tios Sn converging to s), the probability of the Watts event tends to watts(s) 
provided that the distance from c{Dn) to dDn tends to oo. (Oded's SLE con- 
vergence results are similarly general [Lawler, Schramm and Werner (2004), 
Schramm and Sheffield (2005, 2009)].) However, Oded's derivation of Watts' 
formula (like Dubedat's derivation) depends on SLEg convergence, and ex- 
isting SLEg convergence statements [e.g., Camia and Newman (2007)] are 
not quite general enough to imply this. 

Proof of Theorem 7.1. As shown in Section 3, it suffices to prove 
the analogous statement about tripod events {xi,X2), (3:2, X3), and (3:3, X4) 
and the function tripod(s). 

Let £n be a sequence of positive reals tending to zero, and define 7^ to 
be the random interface in D obtained from percolation on £.fi times the 
hexagonal lattice, between the lattice points closest to X2 and X4. Let a„, 
bn, Cn be the tripod set for this interface and the points xi and X3, that is, Cn 
is the first point on 7„ n dDn outside the boundary segment (xi,X3), and 
the interface 7^ up to point c„ last hits the boundary intervals (xi,X2) and 
(x2,a^3) at On and 6^, respectively. From the work of Camia and Newman 
(2007), we can couple the 7„ and 7 in such a way that 7„ — )• 7 almost 
surely in the uniform topology (in which two curves are close if they can 
be parameterized in such a way that they are close at all times). By the 
compactness of dD (and the corresponding compactness — in the topology of 
convergence in law — of the space of measures on dD) it is not hard to see that 
there must be a subsequence of the n values and a coupling of the 7„ with 7 in 
which the entire quadruple (7^1 ^n, ^n, c„) converges almost surely to some 
limit. If we could show further that this limit must be (7,0, &,c) almost 
surely, this would imply the theorem, since uniform topology convergence 
would imply that if 7 hits a before b then 7„, hits a„ before bn for large 
enough n almost surely. However, it is not clear a priori that this limit is 
(7,a, 6, c) almost surely (even though the 7^ converge to 7), since while 7 
touches the boundary at a, b and c, it could be that 7„ comes close to the 
boundary at these points without touching it. 

To obtain a contradiction, let us suppose that there is a uniformly positive 
probability (i.e., bounded away from as n — )• 00) that, say, the limit of 
the On is not a. (The argument for the bn and the c„, is essentially the same.) 
Then there must be an open interval (01,02) of the boundary and an open 
subinterval (/3i,/32) C (ai,a2) of the boundary (with /3i 7^ ai and /32 / 02) 
such that there is a uniformly positive probability that a lies (/3i, /32) but the 
limit of the a„ does not lie in that (01,02)- Now we can expand the Jordan 
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domain D to a larger Jordan domain D that includes a neighborhood of 
(/3i,/52), but where the boundary of D agrees with boundary of D outside of 
(01,02)- Let 7„ denote the discrete interfaces in this expanded domain. We 
can couple the 7„, with the 7„ in such a way that the two agree whenever 7^ 
does not leave D (by using the same percolation to define both). But now we 
have a coupling of the 7^ sequence with the property that there is a positive 
probability that the limit of the 7^ is a path that hits the boundary of 
D\D without entering D\D. This implies that if the jn converge in law, 
they must converge to a random path that with positive probability hits 
the boundary D \ D without entering D \ D. By the Camia-Newman 
theorem, applied to the domain D, the 7„ converge in law to chordal SLEg 
in D, and on the event that SLEg hits d{D \D), it will a.s. enter D\D, 
a contradiction. □ 
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